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Abstract. We analyze the details of mass exchange in the vicinity of the Coulomb barrier for heavy-ion 
collisions involving neutron-rich nuclei using the time-dependent Hartree-Fock (TDHF) theory. We discuss 
the time-dependence of transfer and show that the potential barriers seen by individual single-particle 
states can be considerably different than the effective barrier for the two interacting nuclei having a single 
center-of-mass. For this reason we observe a substantial transfer probability even at energies below the 
effective barrier. 

PACS. 21.60.Jz Nuclear Density Functional Theory and extensions - 24.10.Cn Many-body theory 
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1 introduction 

Heavy-ion fusion reactions are a sensitive probe of the size, 
shape, and structure of atomic nuclei as well as the colli- 
sion dynamics. With the increasing availability of radioac- 
tive ion-beams the study of fusion reactions of neutron- 
rich nuclei are now possible [Tl|^[3]. Other experimental 
frontiers are the synthesis of superheavy nuclei in cold 
and hot fusion reactions ^4u6ii5ii7j|8j , and weakly bound 
light systems [^ [TIIIITTIIT^ . Microscopic descriptions of nu- 
clear fusion may provide us with a better understanding of 
the interplay between the strong. Coulomb, and the weak 
interactions as well as the enhanced static and dynamic 
correlations present in these many-body systems. 

Recently, two aspects of the collision dynamics leading 
to fusion that involve pre-compound neutrons have been 
of interest. Over the last decade a number of fusion studies 
have reported that the average number of neutrons evap- 
orated by the compound nucleus is considerably less than 
what is predicted by statistical fusion evaporation calcu- 
lations '13]. This phenomenon is quite possibly linked to 
the excitation of the pre-compound collective dipole mode, 
which is likely to occur when ions have significantly dif- 
ferent N/Z ratio, and is a reflection of dynamical charge 
equilibration. This was studied in the context of TDHF in 
Refs. [MllTSl lTB] and we have recently observed this phe- 
nomenon in the ^''Ni-l-^^^Sn system [T7]. Similarly, consid- 
erable attention has been given to the influence of neutron 
transfer on fusion cross-sections. Studies suggest that the 
transfer of neutrons with positive Q value strongly en- 
hances the fusion cross-section in comparison to systems 
having negative Q value p!8 l [T9 l [20] . This may explain the 
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fact that lowering of the potential barrier for neutron- 
rich systems does not always lead to higher fusion cross- 
sections. In Ref. [19\ near-barrier fusion of neutron-rich 
nuclei was studied within a channel coupling model for in- 
termediate neutron rearrangement using a semi-empirical 
time-dependent three-body Schrodinger equation. Stud- 
ies showed that for the ''°Ca-|-^^Zr system neutrons were 
transferred in the early stages of the collision from the 
2^5/2 state of ^^Zr to the unoccupied levels of the '^'^Ca 
nucleus. 

It is generally acknowledged that the TDHF method 
provides a useful foundation for a fully microscopic many- 
body theory of low-energy heavy-ion reactions 21j. His- 
torically, fusion in TDHF has been viewed as a final prod- 
uct of two colliding heavy-ions, and the dynamical details 
influencing the formation of the compound system have 
not been carefully dissected in terms of the pre-compound 
properties. Due to the availability of much richer fusion 
data and considerable advances in TDHF codes that make 
no symmetry assumptions and use better effective interac- 
tions, it may now be possible to examine these effects more 
carefully. In TDHF complete fusion proceeds by convert- 
ing the entire relative kinetic energy in the entrance chan- 
nel into internal excitations of a single well-defined com- 
pound nucleus. The dissipation of the relative kinetic en- 
ergy into internal excitations is due to the collisions of the 
nucleons with the "walls" of the self-consistent mean-field 
potential. TDHF studies demonstrate that the randomiza- 
tion of the single-particle motion occurs through repeated 
exchange of nucleons from one nucleus into the other. Con- 
sequently, the equilibration of excitations is very slow, and 
it is sensitive to the details of the evolution of the shape of 
the composite system. This is in contrast to most classical 
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pictures of nuclear fusion, which generally assume near in- 
stantaneous, isotropic equilibration. The relaxation of the 
final compound system is a long-time process occurring 
on a time scale on the order of a few thousand fm/c. In 
contrast, the pre-compound stage corresponds to a time 
scale of a few hundred fm/c. 

In this manuscript we focus on the analysis of trans- 
fer during the early stages of the collision. In particular, 
we confirm the findings of Ref. [19] . We also show that in 
TDHF different single-particle states seem to see differ- 
ent potential barriers in comparison to the generic ion-ion 
barrier. This influences the overall dynamics leading to 
fusion and consequently the effective potential barrier. 



2 Transfer in TDHF 

The TDHF calculations have been carried out using our 
new three-dimensional unrestricted TDHF code [22]. For 
the effective interaction we have used the Skyrme SLy4 
force [23], including all of the time-odd terms. Static 
Hartree-Fock calculations for all the nuclei studied here 
produce spherically symmetric systems. The chosen mesh 
spacing was 1 fm in all three directions, which yield a 
binding energy accuracy of about 50 keV in comparison 
to a spherical Hartree-Fock code. For these calculations 
we have in addition required that the fluctuations in en- 
ergy be as low as 10~'*-10~^, the corresponding accuracy 
in binding energy is about 10^^^. This ensures that the 
tails of the wavefunctions are well converged in the nu- 
merical box. The box size used was 60 fm in the direction 
of the collision axis and 30 fm in the other two directions. 
The initial nuclear separations were 25 fm. 
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Fig. 1. Potential barrier, V{R), for the ^'^O+^'^O system 
obtained from density constrained TDHF calculations (black 
curve). Also shown is the point Coulomb potential (red curve). 



2.1 160+240 system 

As an example of a collision involving one neutron rich 
nucleus we studied the ^^O+^^O system. In order to de- 
termine the potential barrier for the system we have used 
the DC-TDHF method as described in Ref. [21]. In this 
approach the TDHF time-evolution takes place with no 
restrictions. At certain times during the evolution the in- 
stantaneous density is used to perform a static Hartree- 
Fock minimization while holding the neutron and proton 
densities constrained to be the corresponding instanta- 
neous TDHF densities. Some of the effects naturally in- 
cluded in the DC-TDHF calculations are: neck formation, 
particle exchange, internal excitations, and deformation 
effects to all order. The heavy-ion potential was obtained 
by initializing the system at -Ec.m. = 9.5 MeV, which is 
slightly above the barrier shown in Fig. [TJ The peak of 
the barrier is about 8.4 MeV at a nuclear separation of 
9.9 fm. This is lower than the barrier of the ^^O+^^O 
system, which has a height of about 10 MeV. Here and in 
the following, the heavy-ion interaction potential has been 
calculated with a constant mass parameter corresponding 
to the reduced mass of the ions. This is a good approxi- 
mation as long as one is only interested in the value of the 
potential barrier height (as is the case here). For the cal- 
culation of sub-barrier fusion cross sections, however, it is 
essential that coordinate-dependent mass parameters be 
utilized [17] because the cross sections depend sensitively 
on the shape of the potential in the interior region. 

In order to examine the center-of-mass energy depen- 
dence of mass exchange below and above the barrier we 
have initiated TDHF collisions at energies -Ec.m. = 6, 7, 8, 
9, and 9.5 MeV. Interestingly, the head-on (zero impact 
parameter) TDHF collisions for the lowest four energies 
behave as a typical sub-barrier collision: the two ions ap- 
proach a minimum distance with no visible overlap, then 
recoil and move away from each other. This is also true for 
^^c.m. = 9 MeV despite the fact that this energy lies above 
the ion-ion barrier. This suggests that while we can talk 
about an effective ion-ion barrier the individual single- 
particle states may see a barrier somewhat different than 
the effective one. This is in agreement with the findings of 
Ref. [19] , and we shall come back to this point again later 
in the manuscript. Even though we are dealing with sub- 
barrier energies, we observe mass exchange (mainly neu- 
tron) from ^^O to for all these energies. Similarly, we 
observe a small dissipation of the relative kinetic energy, 
ranging from 0.05 MeV for E'en. = 6 MeV to 0.40 MeV 
for Ec.m. = 9 MeV. For comparison, the difference in to- 
tal energy of the system before and after the collision is 
on the order of 0.02 MeV, i.e., the numerical error in to- 
tal energy conservation is less than the dissipated energy, 
even at very low energies. At i?c.m. ~ 9.5 MeV the system 
fuses. 

To gain further insight into the mass exchange we have 
looked at the change in nuclear density along the collision 
axis (z-axis). We define 



p{z,t)= dx dy p{x,y,z,t) , (1) 
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Fig. 2. Partially integrated neutron densities calculated from 
Eq.([l} for the ^'*0 nucleus plotted on a logarithmic scale ver- 
sus the collision axis coordinate z for the ^''0+'^''0 system at 
three energies, -Ec.m. = 7, 8, and 9 MeV. The black-solid curves 
correspond to the initial partial density, the red-dashed curves 
are the same quantity at the distance of closest approach, and 
the blue-solid curves are partial densities long after the recoil. 
Filled spheres near the bottom axis approximately show the 
initial and final location of the two nuclei. 



where p{x,y,z,t) is the instantaneous TDHF density. In 
Fig. [2] we plot p{z,t) for three energies, -Ec.m. = 7, 8, and 
9 MeV, on a logarithmic scale versus the collision axis 
coordinate z. We emphasize that the plotted density con- 
tains only states that correspond to the ^^O nucleus. The 
black-solid curves correspond to the initial partial density, 
the red-dashed curves are the same quantity at the distant 
of closest approach, and the blue-solid curves are partial 
densities long after the recoil. Filled spheres near the bot- 
tom axis denote the approximate initial and final location 
of the two nuclei. The distance of closest approach for en- 
ergies Ec.m. =7, 8, and 9 MeV are 13.2 fm, 11.6 fm, and 
10.1 fm, respectively. Two features are worth noting: first, 
we see a buildup of the density on the ^^O side of the re- 
action plane, and secondly, we note that the buildup gets 
substantially larger with increasing energy. We also ob- 
serve that the recoiled density profile remains very close 
to the initial profile except in the tail region. The cor- 
responding mass transfer AA at sub-barrier energies can 
be found by integrating p{z,t) over the z values on the 
left side of the minimum. At energies Ec.m. = 7, 8, and 
9 MeV we find mass transfer values of 0.024, 0.1464, and 
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Fig. 3. Mass transfer probability defined in Eq. ((Jjl for the 

transfer in the ^®0-f ^'*0 system. 

0.6768, respectively. The mass transfer at E'c.m. — 6 MeV 
is 0.0022. We emphasize that unitarity (or total mass num- 
ber) was conserved with an accuracy of about one part in 
lO'^. In Fig. [3] we plot the energy dependence of the trans- 
fer probability defined as 

where AA is the mass exchange value mentioned above 
and Ain is the mass of the incoming nucleus, in this case 
24. We observe that for sub-barrier energies the curve is 
essentially linear and jumps to one when fusion occurs. 

In order to identify which of the single-particle states 
are predominantly responsible for transfer at sub-barrier 
energies we define the quantity 

Paiz,t) = J dx j dy \^pa{x,y,z,t)f , (3) 

where the ipa^s are single-particle states. In Fig.|4]we plot 
this quantity on a logarithmic scale for the neutron states 
of the ^'^O nucleus at -Ec.m. = 7 MeV, long after the recoil, 
at which time the ion-ion separation is about R = 25 fm. 
Although our single-particle wave functions are calculated 
on a 3-D Cartesian grid and thus do not carry the same 
good quantum numbers as the spherical representation, 
for a well converged spherical nucleus it is possible to cal- 
culate the expectation values of orbital angular momen- 
tum, spin, and parity to identify these states. As expected, 
the contribution to the sub-barrier mass transfer is pri- 
marily coming from the 2s 1/2 neutron state of the 
nucleus. The contribution from the ^5/2 and ^3/2 states is 
an order of magnitude lower than that of the 2si^2 state. 
This again suggests that the barrier seen by the 2si/2 state 
is different than the effective barrier for the entire nucleus. 
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Fig. 4. Neutron single-particle probability densities given by 
Eq. (O for the ^''O nucleus at Ec.m. = 7 MeV at t = Tfinai- 
Filled spheres near the bottom axis approximately show the 
final location of the two nuclei. 
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Fig. 5. Potential barrier, V{R), for the """Ca+^^Zr system 
obtained from density constrained TDHF calculations (black 
curve). Also shown is the point Coulomb potential (red curve). 



2.2 40Ca-|-96Zr system 

In Ref. il9j fusion for neutron-rich systems was stud- 
ied in the vicinity of the Coulomb barrier using a semi- 
empirical time-dependent three-body Schrodinger equa- 
tion. Their studies showed that for the ■*'^Ca-|-^^Zr system 
at £"0.111. — 97 MeV, which is in the vicinity of the Coulomb 
barrier, neutrons were transferred in the early stages of the 
collision {R =11-14 fm) from the 2^5/2 state of ^®Zr to the 
unoccupied levels of the *°Ca nucleus. Here, we shall study 
the same system using the TDHF theory. We have used the 



DC-TDHF method to find the effective potential barrier 
for this system as shown in Fig. [5l The barrier was calcu- 
lated by initializing the TDHF run at £;c.„i. = 97 MeV. 
The barrier peak is about 95 MeV and is located at about 
R= 11.5 fm. For this energy the TDHF collision results in 
fusion. In order to establish the early mass exchange we 
have plotted the nuclear density and the corresponding 
neutron density of the ^^Zr nucleus (which is on the right 
half of the collision plane) in Fig. [6l The neutron density 
was plotted on a logarithmic scale to emphasize the low- 
density contours. The neutron tail seen in the lower frame 
of Fig. [6] is very similar to the neutron density contours 
shown in Fig. 10 of Ref. [T9j. For the same energy we have 
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Fig. 6. Contours of the total density for the dinuclear system 
^"Ca-l-^^Zr at ion-ion separation oi R = 12.3 fm (top frame), 
and contours of the ®®Zr neutron states, plotted on a logarith- 
mic scale at the same separation (bottom frame). 



also calculated the partially integrated neutron density as 
a function of the collision axis coordinate z using Eq.(IT]) 
for the ^^Zr nucleus, shown in Fig. [71 In this case we do not 
have a recoiled final state since the system actually fuses. 
The integrated neutron transfer at the ion-ion separation 
distance of about i? = 12 fm is approximately 0.5. 

To identify which states actually contribute to the 
probability of mass exchange from ^^Zr to ''''Ca we have 
plotted the individual neutron single-particle probabili- 
ties given by Eq. ^ for the ^^Zr nucleus as a function of 
collision coordinate axis z at the same ion-ion separation 
R = 12.3 fm, as shown in Fig.[8l Again, the static Hartree- 
Fock calculation for the ^^Zr nucleus helps us identify 
these states. The states up to = 50 can be enumer- 
ated (by using parity, degeneracy of eigenvalues, etc.) and 
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Fig. 7. Partially integrated neutron density calculated from 
Eq.([l} for the ^^Zr nucleus plotted on a logarithmic scale versus 
the collision axis coordinate z for the '^''Ca+^^Zr system at 
Ec.TCL. ~ 97 MeV. The solid black curve corresponds to the 
initial partial density and the dashed red curve to the same 
quantity at _R = 12.3 fm. 
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Fig. 8. Neutron single-particle probability densities given by 
Eq. © for the ^®Zr nucleus at Ec.m. = 97 MeV at 71 = 12.3 fm. 



exactly match the shell-model states vifith the spin-orbit 
term. The last six states (we do not impose time-reversal 
invariance and thus have one state per nucleon) are al- 
most degenerate in energy and have positive parity values. 
This is consistent with Ref. [19], which finds the neutron 
to be in the 2^5/2 state. The corresponding three states 
(time-reversed pairs) are the red, cyan, and magenta col- 
ored curves drawn thicker than others and pointed out by 
arrows. We observe that while two of the states are the 
largest contributors to the transmission probability, one of 
the states practically makes no contribution. This is true 



despite the fact that all these states are the highest energy 
states and are degenerate in energy. After the second high- 
est 2dt^/2 state at -7.05 MeV (cyan) the next three states 
that have the largest contribution to the transmission 
probability are; one of the two 2^3/2 states at -18.06 MeV, 
the 2pi/2 state at -16.17 MeV, and only one of the five 
1.99/2 states at -12.26 MeV. The single-particle quadrupole 
moment is a good indicator of which of the states, spe- 
cially sub-states belonging to the same quantum number, 
make the largest contribution to the transmission prob- 
ability, since the most stretched states have the largest 
positive quadrupole moments. For example the three sub- 
states of the 2d^/2 state have single-particle quadrupole 
moments of 11.8, 3.7, and —15.6 fm^ corresponding to the 
three states shown in Fig.[8]in descending order. Similarly, 
the 2p3/2 state with a quadrupole moment of 8.0 fm^ is 
the next largest contributor whereas the other sub-state 
with a quadrupole moment of —8.0 fm^ has a transmis- 
sion probability which is orders of magnitude smaller. It 
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Fig. 9. Neutron single-particle probability densities given by 
Eq. Q for the ^"^Zr nucleus at Ec.m. = 91 MeV and after the 
recoiling nuclei are at a separation of i? = 19 fm. 



is interesting to note that some of the sub-states of the 
same quantum number give a much larger contribution, 
and states like I/5/2, which is higher in energy than 2^3/2, 
do not contribute appreciably. This seems to indicate that 
the transmission probability depends not only on the en- 
ergies of the single-particle states, but that it has an addi- 
tional strong dependence on other properties of the states. 
We have also repeated the same study for the '"'Ca-l-^^Zr 
system at Ec.m. = 91 MeV, which is below the effective 
barrier. In this case the nuclei recoil with a closest ap- 
proach distance of 12.6 fm. Again, in order to identify 
which states actually contribute to the probability of mass 
exchange from ^^Zr to ^°Ca we have plotted the individual 
neutron single-particle probabilities given by Eq. ([3]) for 
the ^^Zr nucleus as a function of collision coordinate axis 
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z when the recoiled ions are about i? = 19 fm apart, as 
shown in Fig.O We observe that the single-particle trans- 
mission probabilities after the recoiled nuclei are well sep- 
arated are similar to the findings for the -Ec.m. = 97 MeV 
case, except with reduced probabilities. In addition, the 
observed sub-barrier behavior for the ''°Ca-|-®^Zr system 
is analogous to the ^"^O system. 



3 Conclusions 

We have performed a detailed analysis of mass exchange 
in the vicinity of the Coulomb barrier for systems in- 
volving a neutron-rich nucleus using the TDHF theory. 
Our work was motivated by Ref. [T^ where the same 
phenomenon was studied using a quantum mechanical 
three-body model. For the ''°Ca-|-^^Zr system at i?c.m. — 
97 MeV, which is slightly above the barrier, we can es- 
sentially confirm the results of Ref. ,19J. We confirm that 
at relatively large ion-ion distances the neutron transfer 
probability begins to build up at the location of the re- 
ceiving nucleus. We have analyzed this aspect of neutron 
transmission in TDHF using the lighter ^^0-1-^^0 sys- 
tem. We find that in the vicinity of the Coulomb barrier 
and below the barrier different single-particle states see 
barriers that may differ from the effective ion-ion barrier 
with a fixed center-of-mass. This reflects the fact that the 
single-particle wave functions feel the barrier in the mean 
field potential between the two nuclei, both concerning its 
height and its position relative to their own geometric dis- 
tribution, while the ion-ion barrier is a bulk effect. Due to 
this fact we find appreciable transfer probability for en- 
ergies well below the effective barrier. This confirms that 
the incorporation of neutron transfer effects, as well as 
other effects which depend on the properties of the single- 
particle states, are necessary ingredients in fusion barrier 
calculations. 
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ment of Energy under grant No. DE-FG02-96ER40963 
with Vanderbilt University, and by the German BMBF 
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